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Abstract 

A classical statistical model for the antiferromagnetic (AFM) ordering of 
the Cu-spins in the CuC>2 planes of reduced YBa2Cu306+ x type materials 
is presented. The magnetic phases considered are the experimentally ob- 
served high-temperature AFI phase with ordering vector Qi = (|^0), and 
the low-temperature phases: AFII with Qn = (555) an d intermediate TA 
(Turn Angle) phases TAI, TAII and TAIII with components of both ordering 
vectors. It is shown that the AFII and TA phases result from an effective 
ferromagnetic (FM) type coupling mediated by free spins in the CuO^ basal 
plane. Good agreement with experimental data is obtained for realistic model 
parameters. 
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It is well-established that the C11O2 layers in the high-temperature superconductor ma- 
terials, like La2Cu04 and YBa2Cu306+ x (YBCO), form an antiferrromagetic (AFM) state 
when they are insufficiently hole doped for superconductivity [|I|]. The nearly 2d quantum 
state of these magnetic phases makes them intrinsically interesting, but a major motive for 
studying their properties is clearly the unsettled question, whether the magnetic fluctuations 
contribute to establish the superconducting state or compete with and eventually destroy 
it. 

Despite the many detailed studies of the static and dynamic properties of the magnetic 
phases in YBCO type materials, the origin of the low-temperature AFII magnetic phase has 
not been satisfactorily understood. The AFII phase has been found in the same range of 
oxygen stoichiometries (x < 0.35) as the high temperature AFI phase in apparently pure 
YBCO pm and the homologous rare earth (RE) substituted REBCO However, recent 
experimental studies have shown that it is not present in high-purity single crystals of YBCO 
and NdBCO §, whereas it is found in Al §g, Co and Fe doped YBCO, and 
in YBCO with La doped onto the Ba site [IT|. The early observations of the AFII phase 
in nominally pure YBCO single crystals [^,|3| may therefore result from corrosion of the 
alumina crucibles used for the crystal growth. 

Experimentally it has been shown that the Cu magnetic structures in YBCO materials 
are characterized by a simple AFM ordering of the Cu spins within the Cu02 double layers 
as shown in Fig. la. The Cu spins are lying in the plane but the in-plane easy direction 
has not been established yet. Well below the Neel temperature, (T N = 410 K for x ~ 0), 
the Cu 2+ spins in the Cu02 double layers form a rigid collinear 2c? structure. In the AFI 
phase adjacent double layers {i.e. next nearest Cu02 planes) are AFM coupled via the CuO x . 
basal plane. Thus compared to the chemical unit cell, the magnetic unit cell of the AFI 
structure is doubled along the a and b axes, but not along the c axis giving an ordering vector 
Qi = (||0). In the AFII structure the spins on next nearest Cu02 planes are aligned by an 
effective ferromagnetic (FM) coupling via the CuO x basal plane with a magnetic unit cell 



2 



doubled along the c axis, and hence Qn = (\\\) In an intermediate phase, which 

we shall denote the TA {Turn Angle) phase, components of both vectors, Qi and Qn, are 
observed experimentally. It has been debated whether the TA phase is a mixed AFI-AFII 
phase or results from a continuous transformation between the two spin configurations. In 
this letter we establish a classical statistical model that explains the origin of the interaction 
leading to the AFII phase at low temperatures, and the transition to the AFI phase via 
the TA phase. For pure YBCO there is essentially no free magnetic moment on the Cu 
sites in the CuO x basal plane even for the oxygen content close to the transition into the 
superconducting state, (as ~ 0.3), where significant amount of magnetic Cu 2+ ions should be 
present. This may be explained by the formation of randomly distributed singlet spin-pairs 
of Cu 2+ -0 2 ~-Cu 2+ fragments, which like isolated Cu + ions are non-magnetic, and onwards 
shall be referred to as such. Not only shortest, but also longer chain fragments, like CU3O2, 
should be considered as non-magnetic because the holes may redistribute inside the chains 
and leave the spin state of such a "molecule" as a singlet. It is argued that free Cu spins may 
be formed by direct substitution of magnetic ions like Co || and Fe [|Kj for non-magnetic 
Cu, or indirectly when Al 3+ substitutes for Cu ||[7| or RE 3+ for Ba 2+ PJTT|. In the latter 



case, which may occur for light RE, it is likely that additional oxygen is introduced into 
the basal plane to assure charge neutrality and gives rise to additional oxidation of Cu + to 
Cu 2+ . With RE 3+ substitution for Ba 2+ one free Cu 2+ spin is expected to result, whereas 
Al 3+ may give rise to two Cu 2+ spins, e.g. as shown in Fig. la. These two Cu 2+ spins are 
essentially free because they are linked via non-magnetic Al or the 0-0 orbital overlap. 

The free spins may be polarized and establish an effective FM coupling between the rigid 
spin configurations of two adjacent Cu02 double layers at low temperatures. For sufficiently 
high concentrations of free spins the FM coupling will dominate and form the AFII phase, but 
at higher temperatures it gradually becomes ineffective and the ordering will be controlled 
by the AFM couplings mediated by the non-magnetic Cu ions (cf. Fig. la). Thereby the 
AFII structure is transformed into the TA phase and further into the AFI structure via two 
continuous phase transitions. For low concentrations of free spins the AFII phase does not 



form and the TA phase is the stable ground state. Including a fourfold easy-axis anisotropy 
in addition to a well-established in-plane anisotropy for Cu spins in the bilayers first order 
phase transitions to various TA phases, TAI, TAII and TAII may result. 

The total Hamiltonian for the spin system can be represented via the sum 7i to t = 
J2l=oHk, an d the relevant interaction parameters are defined in Fig. 1. TC is the intra- 
plane Heisenberg spin-exchange and 7i\ is the Heisenberg exchange between the Cu spins in 
the nearest Cu02 planes. H2 is the easy-plane anisotropy and H3 reflects the anisotropy for 
spin rotation in the easy-plane. The indirect Heisenberg exchange interactions between Cu 
spins in the next nearest Cu02 planes via non-magnetic Cu ions, and via free spins in the 
basal CuOa, plane, are given by Hamiltonians H4 and 7^ 5 , respectively. We do not include 
the interactions between spins in the CuO x basal plane, which may be estimated to have 
little influence for low concentrations of free spins. 

In principle the model is applicable for the above mentioned dilute doping mechanisms 
and their combinations. However, in the following we shall focus on the properties of Al- 
doped YBCO. We shall consider any basal plane as consisting of free Cu 2+ spins, of non- 
magnetic Cu ions, and of Al-occupied sites with relative concentration, x, y and 5, that 
satisfy the equality x + y + 5 = 1. Non-magnetic Cu and Al ions are supposed to be 
randomly distributed, while free Cu 2+ spins correlate spatially with Al, e.g. as shown in 
Fig. la. The Cu spins in the j'-th bilayer which form a rigid 2d collinear AFM structure are 
characterized by a classical spin Sj, that is the j-th layer order parameter. The justification 
for this assumption is that the values of the intra-plane and inter-plane coupling constants, 
Jq ~ 170 meV in TL and J\ ~ 10 _1 j7o in Tti, significantly exceed all other contributions 
to Htot HUmi (cf- the experimentally observed temperatures of the AFII and TA ordering: 



^afii < ^ta < 20 K in YBCO). To be specific we ascribe Sj to the lattice site p = (0, 0) of 
the bottom plane (b) of the j-th double layer. Then, evidently, the spin-field in the bottom 
plane is Sf(p) = Sjexp(zQp) and in the top (t) plane sf\p) = —Sf (p). Each spin of 
the set {Sj} governed by the main contribution ?4ot = ^0 + 'Hi is still freely oriented, but 
this degeneracy is removed by smaller interactions, that regulate the relative orientations of 



adjacent double layer spins and their orientations with respect to the easy axes. Because 
the energy scale of TC2, and 7i 4 is much below that of the temperature it is allowed to 
replace them by their "mean" -field expressions in the free energy. Thus Ti,2 restricts the 
spins to lie in the basal plane, and H3 leads to the fourfold easy-axis anisotropy given by: 

JF 3 = _x> 4 cos 40 r -> -2N P V A cos 4 0i ( 1 ) 

r 3 

where <f>j is the angle between Sj and one of the fourfold easy axes, which, in lack of 
information we take to be the x (or equivalently the y) axis. N p is the total number of Cu 
sites in a Cu0 2 plane. Further, 

^4 = J" 2 E sfWflip') - -yN P J2 E S:S j+1 (2) 

3,p' 3 

where the sum runs over the double layers, j, and over the N p y in-plane sites pi which are 
bridged by non-magnetic Cu ions in the basal plane. 

The Hamiltonian 7i 5 realizes coupling via free Cu spins: 

«s = ^E<^(P")(S?V) + Sfi(p")) (3) 

3,P" 

where now the summation is over j and the xN p sites p" which are bridged by free Cu spins, 
cfj + i- The corresponding free energy obtained from a classical statistical average over the 
angular degrees of freedom of <x- + i is: 



jr 5 = - x N p k B T\n 
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which behaves as —xN p \^F 3 \ a \ Sj — Sj + i\ when T — > and goes to zero for T — > 00. Thus we 
arrive at the effective free energy, jF eff = ^=3 -^i? which takes the form of an effective Id spin- 
system {Sj}. Then our goal is to minimize jF eff . The polarization type of coupling leading 
to JF 5 tends to align the order parameters Sj and S J+ i in antiparallel at low temperatures, 
whereas it becomes ineffective at high temperatures, where the parallel alignment of Sj and 
Sj + i will prevail, due to JF 4 if J 2 > 0. At intermediate temperatures the competition between 
AFI (parallel) and AFII (antiparallel) alignments of the spins in adjacent bilayers will lead 
to the TA phases. To show this we note, that the in-plane angles may be presented by 
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(j>j = + (-l)ty (5) 

where ip is the turn angle and ipo = n\ . In i/; oddn (TAII phase) as well as evenn (TAI and 
TAIII phases) are necessary to establish all the minima of JF 3 . Fig. lb illustrates the various 
possibilities, including the pure states, AFI and AFII, corresponding to ip equal to and 
f, respectively. Note, that Eq. (j5|) defines the spin rotation angles as strictly alternating: 
A(pj = (f)j + i — <pj = (— l)^ 1 ■ 2ip. For £> 4 = 0, this is not justified because A0j may take 
the form of degenerate configurations Rj ■ 2ip with Rj being a random set of ±1. As a 
result the spin configuration becomes disordered along z. A similar disordering occurs in 
the TAII phase for ip = |, which means A(f>j — Rj • | and is just the case shown in Fig. la. 
In principle, longer range interactions prevent the above mentioned disordering to happen. 
■ip = I or y would also be candidates for disordering but because they appear at first order 
transitions between TAII and TAI or TAIII phases, this kind of a disordering is inaccessible. 
Thus we arrive at the normalized free energy functionals (f^ 1 = J r e ff/(2N p \J' 3 \Sa)): 

f±(z) = yjz* ± g(, 2 - z*) T f - xt In (6) 

where the normalized parameters are: 

. J2S 8E> 4 k B T 

and / + and / _ refer to the turn angles centered around the easy (TAI and TAIII phases) 
and hard (TAII phase) directions, respectively. 

For the concentration of free spin, x, we shall assume that each Al 3+ ion generates two 
free Cu spins, i.e. x = 25 (cf. Fig. la) and thereby: y = 1 — |x. 

f ± {z) have been minimized with respect to z as function of parameters, t (reduced 
temperature) and x, for a few values of the normalized interaction parameters j and q. A 
realistic estimate for Ji is ~ 0.1 K. Much less is known about 3% and P4. We shall leave T>± 
as a free parameter. 3% cannot be too small, because like in the planes this spin interaction 
is mediated via oxygen ions. However, hybridization via apical oxygens is not as effective as 
hybridizations within the planes or chains because the involved orbitals are ortogonal. Since 



the AFI - AFII transition temperature is of order a few Kelvin we estimate as 10 K and 
thereby: j = 0.01 and T = 5t K. 

Fig. 2 shows the corresponding phase diagram for three values of the anisotropy parame- 
ter q. Full lines represent continuous transitions, dotted lines are first order transitions, and 
dots are tricritical points. The dashed lines are the disorder line at z = |. Although the 
anisotropy q does not change the phase diagram boundaries significantly, it influences the 
nature of the phase transitions even when q -C j . The AFI-TAI transition is always continu- 
ous, and if q = 0, the AFI-AFII transition always takes place via two continuous transitions 
through the intermediate TA phase. For low concentrations of free spins (x<2jy), the pure 
AFII state is inaccessible even at t — > (see Fig. 2). Increasing j to 0.1 does not change 
the morphology of the phase diagram but for obvious reasons higher concentrations of free 
spins are required for the characteristic features to occur. One can subdivide the order 
parameter in the TA phases into AFI and AFII parts as ttiafi = §|Sj + ^i+il = Sy 1 — Zg 
and ?7i afii = -| | Sj — Sj + i\ = Sz g . z g is the z value at which the global minimun of the free 
energy is achieved from f^ 1 in Eqs.(|6|). The inserts in Fig. 2 show examples of the tempera- 
ture variation of the global order parameter, z g , as composed from the optimal z . Without 
anisotropy (q = 0) continuous transitions are observed, but for q = 0.001 and x = 0.02 a 
continuous transition from AFI (z g = 0) to TAI (z g = z + ) is followed by the first order 
transition to TAII (z g = z~) and then to AFII (z g = 1). 

A direct comparison with experimental data requires a detailed knowledge of the number 
of free spins generated for specific Al doping levels. It is well-established that the sample 
treatment is very essential for the Al configurations in the basal plane, and quite different 
concentrations of free spins may develop for the same Al stoichiometry [[14]]. However, 
comparing the model results displayed in Fig. 2 with experimental data reported for Al 
doped YBCO ||[?| we find agreement with the following observations: 

(1) the sequental transitions AFI-TA-AFII occur below 20 K (i.e. below t ~ 4); 

(2) the AFII and TA phases are not present in nominally pure materials (i.e. with x = 0); 

(3) the transition into the AFII phase is not always complete and the TA phase remains 
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stable at low temperatures and small x; 

(4) there is essentially no ordered magnetic moments on the Cu sites in the basal plane; 

(5) the intensities of magnetic Bragg peaks at Qi and Qn are in good agreement with 
the square of the order parameters ttiafi and 7?iafii, (c/. insert to Fig. 2a). 

Since apparently no first order transitions or jumps in the m^-pi and ttiafii order param- 
eters are observed experimentally, the four-fold easy axis anisotropy is very small. 

When large amount of free spins are introduced in the basal plane, as may be the case 
when YBCO is heavily doped with Co || or Fe l lOfl , they interact and it is likely that 
they form a 2d AFM structure. In this case the AFII phase is expected to be stable even 
at higher temperatures as has been observed experimentally in Fe doped YBCO, where 
complete suppression of the AFI phase may result |10[ . 

A classical treatment of a spin-| system may appear inappropriate. However, it should 
be recalled, that we are mainly giving a qualitative mechanism for the origin of the AFI , 
AFII and TA phases. A quantum treatment will be pursued. 

In conclusion we have shown by a classical statistical model that the AFII and TA phases, 
observed experimentally in Al, Co and Fe doped YBCO and in materials with rare earth 
(RE) ions on the Ba site, result from polarization of free spins created in the basal plane of 
the structure. Thus, we propose that in nominally pure materials with no RE on the Ba site 
the AFII phase should not be present. Using realistic values for the interaction parameters 
a semi-quantitative agreement with the experimental observations is established. Based on 
the experimental fact that the magnetic transitions appear to be continuous and long range 
order is etablished along the z axis we conclude that the in-plane anisotropy parameter, T>^ 
is finite but very small (< 10 5 1 1 , cf. Fig. 2b). More model work and experimental studies 
on carefully prepared materials are needed and in progress for more detailed comparisons. 
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FIG. 1 Magnetic structures observed in YBCO. (a) Interactions Jq and J\ are responsible 
for the rigid double layer magnetic structures; Sj is the bilayer order parameter; J2 the coupling 
via non-magnetic Cu + that favors AFI structure, while J$ via the free Cu 2+ spins tends to AFII 
ordering. The formation of Cu 2+ spins through Al substitution is also shown. The mutual spin 
arrangement in the bilayers is related to the TAII phase, (b) displays spin configurations and 
turn angles (shaded areas) tp in adjacent layers j and j + 1; magnetic phases AFI, AFII as well as 
various TA phases may result from the model studies. are related to hard and easy axis spin 
anisotropies. 
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FIG. 2 Magnetic phase diagrams for Al doped YBCO as function of the free spin concentration 
x and the reduced temperature, t, for the normalized interaction constant j = 0.01 and different 
values of the anisotropy q. The diagram for q = is representative of a finite but infinitely small 
anisotropy (see the text). The dashed lines are the disorder line as explained in the text. Inserts 
show the order parameters for the selected x values (vertical dashed lines). 
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